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Abstract. We report the detection of massive star forma- 
tion along a bar in the peculiar starburst galaxy Mkn 439. 
We present optical B, R and Ha+[NII] emission line im- 
ages as well as H band images to show that the signature 
of the bar becomes progressively weak at longer wave- 
lengths. Moreover, this bar is misaligned with the main 
body of the galaxy. The peak Ha emission does not co- 
incide with the bluest regions seen in the colour maps. 
We infer that the starburst is young since the stars in the 
burst have not started influencing the light in the near 
infrared. There are indications of dust in the inner regions 
of this galaxy. 
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tailed surface photometric studies of this galaxy have been 



reported. Usui, Saito & Tomita (1998) report the detec 
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1. Introduction 

Mkn 439 (NGC 4369, UGC 7489, IRAS 12221+3939) is 
a nearby early type starburst galaxy (z=(h0035). It has 
been classified as a starburst by Balzano ( 1985 ) based on 
the equivalent width of Ha and als o belo ngs to the IRAS 
Bright Galaxy Sample (Soifer et al. 1987 ). On the basis of 
multiaperture near i nfrare d photometry and optical spec- 
troscopy, Devereux ( 1989 ) describes t his ga laxy as a M82 
type starburst galaxy. Rudnick & Rix ( 1998 ) report an az- 
imuthal asymmetry in the stellar mass distribution of Mkn 
439 based on the R band surface brightness. The peculiar 
morphology of Mkn 439 attracted our attention during the 
course of an optical imaging study of a sample of starburst 
galaxies derived from the Markarian lists (Chitre 1999| ). 
The galaxy image was nearly circular and appeared fea- 
tureless in long exposure images. The outer isophotes were 
smooth and nearly circular in B and R bands. However, 
the isophotal contours show highly complex features in 
the inner parts. Moreover, the strength of these features 
is wavelength dependent. Wiklind & Henkel (1989) report 
the detection of a molecular bar in the central region of 
this galaxy based on the results of CO mapping. No de- 



tion of two regions bright in Ha that are displaced from 
the nucleus and faint emission from the nucleus. However, 
their data were obtained at a seeing of 5". In order to 
study the spatial distribution of various stellar popula- 
tions in Mkn 439, we imaged this galaxy in B, R, Ha and 
H bands. The B and R band continuum trace the interme- 
diate age populations while Ha traces the young, massive 
stellar populations. The infrared continuum of galaxies is 
dominated by evolved stellar populations. Hence, the H 
band and the line emission images can be used alongwith 
the optical continuum images to separate the young and 
old stellar populations spatially. 



2. Observations and data reduction 

2.1. Optical (B,R) and Ha imaging 

The B, R and Ha images were obtained under photomet- 
ric conditions from the 1.2m telescope at Gurushikhar, 
Mt. Abu. The images were taken at the Cassegrain fo- 
cus employing a thinned back illuminated Tektronix IK 
x IK CCD. Binning of 2 x 2 was employed before record- 
ing the images to increase the signal-to-noise ratio of the 
measurements and keeping in mind the data storage re- 
quirements. The final resolution was 0.634"/pixel which is 
sufficient to sample the point spread function (PSF) ap- 
propriately. Typical seeing (full width at half maximum 
(FWHM) of the stellar images) was ~ 1.8" for the im- 
ages. For the Ha images a filter having FWHM of 80 Awas 
used. Another off-band filter of the same FWHM was used 
to measure the galactic red continuum. About 3-4 expo- 
sures were taken in each of the photometric bands. The 
total exposure times were 510 sec, 360 sec and 1600 sec in 
B, R and Ha respectively. Standard stars from Landolt 



(1992) were observed to calibrate the broad band data. 
Twilight flats were taken and median filtered to construct 
the master flats. The data was reduced using IRAF [] on 
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1 IRAF is distributed by National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities 
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Table 1. Global properties of Mkn 439 



Parameter 



Value 



a (2000) 


12 ft 22 m 0.8 s .4 


5(2000) 


39° 39 41 


RC3 type 


RSAT1 


UGC type 


SO/Sa 


b Adopted distance 


18 Mpc 


a R ° 

JDrp 


12.27 


a (f/-B) T 


-0.02 


& {B-V)t 


0.65 


Lp/ij 


4xl0 9 L Q 


a RC3 




b Deutsch & Willner (1987) 



the IBM-6000 RISC at PRL, Ahmedabad. A detai led re - 
duction procedure can be found in Chitre & Joshi ( 1999[ ). 



2.2. H band images 

The H band images were recorded with a 256x256 NIC- 
MO S array at the 1.2 m Gurushikhar telescope. An expo- 
sure of 30 sec ensured that the background and the galaxy 
signal were in the linea r portion of the detector response 
curve (Joshi et al. 1999). Observations were made by alter- 
nating between the galaxy and positions 4'-5'to the north 
and south till a total integration time of 600 seconds on 
the galaxy was achieved. Several dark frames having the 
same time sequences as that of galaxy or sky were taken 
and median filtered master dark frames were constructed. 
The median filtered master sky frames were constructed 
using several sky frames with integration times equal to 
those given for the galaxy. All the source frames were cor- 
rected for the sky background by subtracting the master 
sky frame from the source frames. As the program galaxy 
does not occupy the whole detector array, the residual sky 
was determined from the image corners and the images 
were then corrected for residual sky. The dark subtracted 
sky frame was used to construct the master flat. The sky 
corrected galaxy frames were corrected for fiat field re- 
sponse of the detector by dividing the galaxy frames by 
the master flat. 

Finally, the galaxy images were aligned by finding the cen- 
ter of the galaxy nucleus using the IMCNTR task in IRAF 
and co-added to improve the S/N ratio. The plate scale 
was selected to be 0".5 per pixel. Faint standard stars 
from the UKIRT lists were observed for calibration. 



3. Morphology of Mkn 439 

Fig. [y illustrates the isophotal contours of the inner 25" 
of Mkn 439 in B, R, continuum subtracted Ha and H 



Inc. (AURA) under cooperative agreement with the National 
Science Foundation, USA. 



band. A comparison of the various panels in Fig. Vu shows 
that the morphological structures vary at different wave- 
lengths. The morphology of Mkn 439 in the B band is 
characterized by smooth outer isophotes and a very com- 
plex light distribution in the inner region. The central 
region is elliptical and is elongated in the NS direction. 
Faint indications of a spiral arm in the NE direction is 
seen in the isophotal maps in B and R. The contour maps 
show two projections - one along the NW and the other 
along the SE from the nuclear region. These projections 
are most prominent in the B continuum, getting progres- 
sively fainter at longer wavelengths and nearly disappear- 
ing in H . The B band image shows another condensation 
to the SW of the nuclear region. Similar to the projections, 
this feature also gets progressively fainter at longer wave- 
length. The H band image shows smoother isophotes. The 
signature of the projections is absent at this wavelength. 
As seen in the R band, the outer isophotes are nearly cir- 
cular. However, unlike other optical bands, there are no 
spurs or bar-like features apparent in the H band image. 
The continuum subtracted Ha image shows an elongated 
bar-like strucure corresponding to the projections seen in 
the contour maps. Ha emission is seen along the bar in 
the form of clumps. Emission is most intense at the ends 
of the bar, though it is found to extend throughout the 
body of the galaxy. Emission from the nucleus is much 
fainter as compared to that from the clumps in the bar 
ends. Ha emission is maximum in Spot 1. The clump of 
Ha emission seen to the E of the extended emission has 
no counterpart in the continuum colour map. The bright 
blobs of emission in Ha have no counterparts in the H 
band. This indicates that the HII regions are young and 
have not yet evolved enough to form a considerable num- 
ber of red giants and supergiants to start influencing the 
light in the H band. It is also seen that the latest episode 
of star formation is misaligned with the isophotal con- 
tours of the near infrared continuum. The (B-H) colour 
map (Fig. |2J) was constructed by scaling the images, ro- 
tating and aligning them. It shows interesting features. A 
bar-like structure made up of blue clumps is seen in the 
central part of the galaxy. A spiral arm starts from the 
nuclear region and curves towards the eastern side. A dis- 
tinct blue clump is present at either end of the bar marked 
as Spot 1 and Spot 2 in Fig. g. These correspond to the 
ends of the two projections seen in the isophotal contours 
in B. Another blue region (Spot 3) is seen about 8" to 
the south of Spot 1. The clump of Ha emission seen to 
the E of the extended emission has no counterpart in the 
continuum colour map. The (B-R) and (B-H) colours of 
these regions are listed in Table 2. The isophotal contours 
of Mkn 439 appear different in the optical, near infrared 
and the line emission indicating the spatial separation of 
the distribution of these various populations. The gaseous 
component in the galaxy appears to be under the influ- 
ence of a potential which has distributed it in the form 
of a gaseous bar. Compression of the gas in the bar has 
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Fig. 1. Top panel: B (left) R (right) band contours. Bottom panel: H (left) and continuum subtracted Ha+[NII) 
(right) isophotal contours of Mkn 439. North is at the top and East is to the left. The isophotal contours are plotted 
on a logarithmic scale at intervals of 0™\3 with the peak contours being at 17 m .5, 17 m .O and 14'™. 7 for B, R and H 
respectively. The Ha+[NII] contours are plotted such that the lowest contour corresponds to 2a of the background. 



led to the formation of young, massive stars which are 
seen as clumpy HII regions along the bar. We infer that 
the latest dynamical episode experienced by the galaxy 
has given rise to the formation of young, massive stars 
along the bar as a result of the response of the gas to the 
perturbing potential. A comparison of the Ha contours 
in Fig. H and Fig. reveals that no HII regions are seen 
in the blue spiral arm-like feature emerging from the nu- 
cleus indicating that the blue spiral arm is made up of 
an intermediate age stellar population. Wiklind & Hcnkcl 



(1989) report the detection of a molecular bar in the cen- 
tral region in this galaxy based on CO mapping. They ob- 



served Mkn 439 in both the J=l-0 and J=2-l line of CO 
and found that the ratio of the J=2-l to the J=l-0 inten- 
sity varies with position and inferred that this was due to 
changing physical conditions in the molecular cloud pop- 
ulation. The contour nraps of these two transitions can be 



found in Wiklind (1990). Many galaxies with weak stellar 
bars have been found to contain pronounced bar-like gas 
distributions similar to the one found in Mkn 439. For ex- 
ample, the center of the nearby Scd galaxy IC 342 harbors 
a bar-like molecular gas structure and a modest nuclear 



starburst (Lo et al. 1984; Ishizuki et al. 199C). Other ex- 
amples of galaxies having a molecular bar at their centers 
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Fig. 2. B-H colour map of Mkn 439 
Table 2. Colours of clumpy regions 



region B-R B-H 

nucleus 04 21S 

spot 1 0.7 2.9 

spot 2 0.6 2.6 

spot 3 0.8 3.2 



are NGC 253 and M83. Simulations by Combes ( J1994J ) de- 
scribe the formation of a gas bar which is phase shifted 
from the stellar component in the innermost regions of a 
galaxy due to the existence of perpendicular orbits. How- 
ever, her models describe the situation for the nuclear bars 
and in the innermost lkpc region. An alternative explana- 
tion could be that two unequal mass spira ls have merged 
to form the SO galaxy. Bekki Kenji ( 1998 ) suggests that 
SO galaxies are formed by the merging of spirals and when 
two spirals are of unequal mass, the SO galaxy thus formed 
has an outer diffuse stellar envelope or a diffuse disk like 
component and a central thin stellar bar composed mainly 
of new stars. 

4. Isophotal analysis 

In order to provide a quantitative description of the mor- 
phological aspects at various wavelengths, we explored 
Mkn 439 using ellipse fitting techniques. The procedure 
consists of fitting elliptical isophotes to the galaxy im- 
ages and deriving 1-dimensional azimuthally averaged ra- 
dial profiles for the surface brightness, ellipticity and the 
position angle based on the algorithm given by Jedrejew- 



ski (1987). This technique has been used successfully in 



studying various structures in galaxies like bars, rings, 
shells, etc. a nd in searching for d ust in them (Bender 
& Mollenhoff [L987t Wozniak et al. |l995| and Jungweirt, 




5 10 15 20 25 



0.4 -f 



>S^>a*,** 



5 10 15 20 25 



150 
100 
50 



D D CT 



5 10 15 20 25 



Fig. 3. The surface brightness, ellipticity (e) and position 
angle (P.A.) profiles of Mkn 439. The fitted disk in the H 
band is shown by the line. 
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Fig. 4. Radial distribution of B-R and B-H colours 



Combes & Axon 1997). Multiband isophotal analysis can 



also be used to indicate whether the reddening seen in 
colour maps is due to a redder stellar population or due 
to the presence of dust (Prieto et al. |1992a| , |1992b| ). The 
surface brightness distribution and the variation of the 
position angle and ellipticity of the isophotes in each fil- 



A.Chitre: Star formation in Mkn 439 



ter (Fig. ||) were obtained by fitting ellipses to the images 
in each filter using the ISOPHOTE package within STS- 
DASn. The detailed fitting procedure used is outlined in 
Chitrc ( |1999| ) . The radial distribution of the colour indices 
(Fig. 0) was derived from the surface brightness profiles. 
Fitting isophotes to the images reveals changing elliptic- 
ity and position angle throughout the body of the galaxy 
(refer Fig 0). The luminosity profile is smooth except for 
small features at 5" and 10" in the optical bands. An in- 
spection of Fig. [| shows that the galaxy is bluest near 
the center and gets redder outwards. The ellipticity of the 
elliptical feature is maximum at the center and goes on 
decreasing outwards unlike a bar in which ellipticity in- 
creases accompanied by a constant position angle. The 
ellipticity profile shows a double peaked structure in the 
inner region. The first peak is seen between 2"- 3" and the 
second peak at 5". The ellipticity of the first peak is wave- 
length dependent, the isophotes at shorter wavelengths 
being rounder. The colour map also shows a small local 
redder region between 3" and 4". The surface brightness 
profiles also show a small dip in the intensity at shorter 
wavelengths at 4". All these features indicate the presence 
of dust in the inner 4" of this galaxy, van den Bergh & 
Pierce ( |1990D do not find any trace of dust in Mkn 439 
from a direct inspection of the B band images on a CCD 
frame. However, ellipse fitting analysis has been success- 
fully employed in the present study to infer the presence 
of dust in the inner regions of this galaxy based on multi- 
band observations. The other peak occurs at 5" which 
corresponds to the brightest region seen in Ha. The depth 
of the dip between the two peak reduces at longer wave- 
lengths. The first peak and the dip is probably due to dust 
while the second one corresponds to the blue region at the 
end of the bar. Both these factors, namely dust and star 
forming regions contribute the maximum at shorter wave- 
lengths. At longer wavelengths, the effects of both dust 
and the star forming regions are reduced hence we see the 
underlying old stellar population. As a result the depth of 
the dip is reduced at longer wavelengths. Beyond 5", the 
ellipticity starts dropping and reaches a value (^0.05) at 
15" and remains at a low value beyond that in all filters. 
Between 5" and 15", the isophotes at shorter wavelengths 
are rounder than the corresponding isophotes at longer 
wavelengths indicating the presence of dust in this region 
of Mkn 439. The position angle is nearly constant in the 
inner 10". The luminosity profiles show an inner steeply 
rising part and an outer exponential disk. We derived the 
scale lengths of Mkn 439 in each of the filter bands. This 
was done by marking the disk and fitting an exponential 
to the surface brightness profile in this region. The range 
of fit was taken to be from 18" to the region where the sig- 
nal falls to 2ct of the background. The fit to the H band is 
shown in FigM The scale lengths derived were 0.97±0.14 



2 The Space Telescope Science Data Analysis System STS- 
DAS is distributed by the Space Telescope Science Institute. 



kpc in B, 0.84±0.02 kpc in R and 0.61±0.03 kpc in H 
band. 

5. Conclusions 

1. Mkn 439 is a peculiar galaxy made up of three distinct 
components: an elliptical structure in the inner regions, 
a smooth outer envelope in which this structure is em- 
bedded and a bar. We detect massive star formation 
along the bar in Mkn 439. This bar is misaligned with 
the main body of the galaxy. 

2. The signature of the bar gets progressively fainter at 
longer wavelengths. 

3. The stars in the bar are young and have not yet started 
influencing the light in the near infrared region. This 
indicates that the galaxy has undergone some pertur- 
bation which trigerred the bar formation and the star- 
burst along the bar in recent times. 

4. There are indications for the presence of dust in the 
inner 15" of the galaxy. 
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